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ABSTRACT 

Modeled dissipation rate transport equations are often derived by invoking various hy- 
potheses to close correlations in the corresponding exact equations. D. C. Leslie suggested 
that these models might be derived instead from Kraichnan’s wavenumber space integrals 
for inertial range transport power. This suggestion is applied to the destruction terms in 
the dissipation rate equations for incompressible turbulence, buoyant turbulence, rotating 
incompressible turbulence, and rotating buoyant turbulence. Model constants like C s 2 are 
expressed as integrals; convergence of these integrals implies the absence of Reynolds num- 
ber dependence in the corresponding destruction term. The dependence of C e 2 on rotation 
rate emerges naturally; sensitization of the modeled dissipation rate equation to rotation 
is not required. A buoyancy related effect which is absent in the exact transport equation 
for temperature variance dissipation, but which sometimes improves computational pre- 
dictions, also arises naturally. Both the presence of this effect and the appropriate time 
scale in the modeled transport equation depend on whether Bolgiano or Kolmogorov in- 
ertial range scaling applies. A simple application of these methods leads to a preliminary 
dissipation rate equation for rotating buoyant turbulence. 


* This research was supported by the National Aeronautics and Space Administration 
under NASA Contract No. NAS1-19480 while the authors were in residence at the Institute 
for Computer Applications in Science and Engineering (ICASE), NASA Langley Research 
Center, Hampton, VA. 



I. Introduction 


Despite widespread agreement in the turbulence modeling community about the basic 
structure of the dissipation rate transport equation, questions remain about its theoretical 
grounding. The central problem remains to justify the absence of 0(-Re?/ 2 ) terms in this 
equation, the possible existence of which is suggested by elementary arguments. 1 More 
basic difficulties arise when there are coupled fluctuating fields as in buoyant turbulence, 
or imposed time scales as in rotating turbulence. In these cases, even the form of the 
appropriate transport equation is doubtful. Thus, it has been suggested 2 that the transport 
equation for dissipation of temperature variance eg should contain a term related to the 
buoyant production of turbulence energy, although this requirement does not follow from 
the exact transport equation for eg. Similarly, although rotation is well known to have 
effects which must be reflected in the dissipation rate transport equation, such effects do 
not occur explicitly in the exact equation for e. These effects are therefore sometimes 
introduced by sensitizing 3 the transport equation to rotation in various ways. 

Attempts to derive the e transport equation from analytical theories have been only 
partially successful. In a preliminary discussion of two equation modeling based on the 
direct interaction approximation 4 (DIA), Leslie 5 observed that whereas the transport equa- 
tion for turbulence kinetic energy followed in principle by integrating the DIA correlation 
equation over all wavevectors, the analogous treatment of the response equation, which 
should lead to a transport equation for a time scale, does not suggest itself as readily. 
Leslie asserted that dimensional analysis alone fixes the functional form of the e equation, 
but this statement neglects the possibility of Rex dependence. A subsequent analytical 
investigation of the dissipation rate transport equation by Yoshizawa 6 led to 0{Re\l 2 ) 
terms and to dependence on the energy spectrum in the dissipation range. A later analysis 
by Yakhot and Orszag 7 was criticized by Smith and Reynolds 8 and then reformulated by 
Smith and Yakhot. 9 This work identified the absence of 0(Re ^ 2 ) terms with a cancellation 
of divergences in the renormalization group formulation of the dissipation rate equation. 
This point of view must be contrasted with the usual one in the modeling literature, namely 
that 0(Re-x 2 ) terms are absent because certain production and dissipation terms in the 
exact transport equation for dissipation rate cancel. The failure of this cancellation, un- 
balanced vortex stretching, has been analyzed by Speziale and Bernard, 10 and of course 
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leads to explicit Rex dependent terms. 

These investigations of the e transport equation follow the usual approach in tur- 
bulence modeling of closing correlations in the exact equation for i —< i/Vu • Vu >. 
A different approach is taken by Hamba, 11 namely that a transport equation for eddy 
viscosity is obtained by Yoshizawa’s two-scale direct interaction approximation 12 ’ 13 and 
the transport equation for e is then inferred from this equation and a K equation. The 
present analysis is based instead on Leslie’s suggestion 5 that the £ transport equation be 
derived from Kraichnan’s integrals 14 for inertial range transport power. Leslie’s original 
discussion was general enough to include effects of production and turbulent diffusion; the 
present analysis is restricted to issues accessible to a homogeneous theory. Applied to 
incompressible turbulence, this procedure leads to an integral for C e 2 ; the convergence of 
this integral at high wavenumbers is equivalent to the absence of Rex dependence of the 
destruction term. A modeled dissipation rate equation for rotating turbulence is derived 
by substituting rotation- dependent field descriptors in Kraichnan’s integrals. To illustrate 
the general method, the phenomenological theory of rotating turbulence of Zhou 15 is ap- 
plied. The dependence of C S 2 on rotation emerges naturally, and even provides a result 
very similar to models proposed by Speziale. 16 Buoyant turbulence is investigated by writ- 
ing the analogous integral for eg. The analysis leads naturally to the buoyancy dependent 
term suggested in Ref. 2. Straighforward application of the arguments of Ref. 15 to 
buoyant turbulence leads to a preliminary dissipation rate equation for rotating buoyant 
turbulence. The present approach to the dissipation rate equation, like that of Hamba and 
Yoshizawa, offers an alternative to the standard derivations. Application of this approach 
to production and diffusion terms will be reported later. 

II. The dissipation rate equation in incompressible turbulence 

It is generally accepted that inertial range spectral transport power must be balanced 
by viscous dissipation; accordingly, equating dissipation to the inertial range transfer in- 
tegral of Kraichnan, 14 

e = % -[I~{k') - I + (k')]P imn ( k) < u m (p)u n (q)ui(-k) > (1) 

where 

Pimn{ k) = k m Pi n (li) -f- &n-Pim(k) 

Pij( k) = Sij — kikj / k 2 
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and the integration operators are defined by 

r 2k' 


I+(k') = / dk / 
Jk' Jp 


dp dq 


p,q<k' ; k=p+q 


■(*') = f / 

Jo Jp 


dp dq 


>p,q>k' ; k=p+q 

Time arguments in the Fourier amplitudes Ui are understood. Differentiate Eq. (1) with 
respect to time, substitute the equations of motion for the u terms, and apply the quasi- 
normal hypothesis to obtain 


£ — -[I {k ) I + (k )]Fj mn (k) {<! Ujn(p)^n(q)^i( k) > + < u w (p)?i n (q)'Uj( k) !> 

+ < u m (p)u re (q)ui(-k) >} 

= \[r(k') - I + (k')]{2P imn (k)P mra (p)Q ns (q)Q ir (k) 

Fjmn(k)i :> 1 i rs (k)Q ns (q)Q mr (p) } (2) 


In turbulence with an isotropic inertial range, the correlation tensor takes the form 


Qij{ k) = Q{k)Pij{ k) 


(3) 


where the correlation function Q is related to the energy spectrum by 


Airk 2 Q[k) = E(k) 


(4) 


The integrals in Eq. (2) will be evaluated by substituting the Kolmogorov spectrum 


E(k) = C K £ 2/3 k~ 5 / 3 


(5) 


The absence of Rex dependence will follow from the convergence of these integrals at high 
wavenumbers: if the integrals diverged at either high or low wavenumbers, it would be 
necessary to substitute instead a truncated Kolmogorov spectrum such as 


E ( k)= fC K eWk-s/> iik 0 <k<k d 

\ 0 if k < k 0 or k > kd 


( 6 ) 


and investigate the dependence on the cutoffs ko and k&.. Dependence on kd- due to 
divergence at high wavenumbers, would imply Reynolds number dependence through the 
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relation Rex = (kd/k 0 ) 4 ^ 3 . The convergence of the integrals Eq. (2) would prove the 
locality of e, namely that e is a property of the inertial range alone. 

Substitute the isotropic forms Eqs. (3)-(5) in Eq. (2). The convergence of each term 
will be investigated separately. 

There are two I~ terms: 

Tf = /-p imn (k)P mrs (p)P ns (q)P ir (k)Q( ? )Q(fc) (7) 

^2 = ^ Pim»(k)Pfr#(k)P) ns(_ f l)Pmr(jp')Q{p)Q{q) (8) 

If p — > oo with k fixed, the triangle condition k = p + q implies that q—>p. First consider 
the integral defined in Eq. (7). Substitute the leading order result q = p; ZJ~ vanishes 
to this order because it is of odd order in p. Expand the q dependent terms to next order 
in k /p. Then Zf is of the order 

pk' poo pk* poo 

Z^ ~ / E(k ) dk / dp E(p)kp(k/p) = / dk k 2 E(k) / dp E(p) 

Jo Jk’ Jo J k’ 

- .(: fc ') 4 / 3 ( fc , r 2/8 = ( fc ') 2/3 ( 9 ) 

For the term I 2 of Eq. (8), the scaling is found directly to be 

pk f poo 

I 2 ~ / k 4 dk / dp p 2 p -22 / 3 ~ (fc , ) 5 (fc , ) -13//3 ~ (A:') 2 / 3 (10) 

Jo Jk' 

To investigate the convergence of the I + terms, note first that since the wavevectors 
are all bounded from above in these integrals, there is no question of high wavenumber 
divergence, so consider their convergence at low wavenumbers. It will be convenient to 
rewrite the I + terms symmetrically as 

Z + = I + Pimn(k){Pmrs(p)Pns{<l)Pir(k)Q(q)Q(k) + P mr s(q)Pns(p)Pir(k)Q(p)Q(k) 

- P ira (k)P ns (q)P m r(p)Q(?)Q(p)} (11) 

It suffices to check convergence in the limit q -+ 0. The triangle condition implies that 
also p — > k. Substitute the lowest order result p = k in first and third terms in Eq. (11) 
to obtain the integrand 

Pimn (k)P res (q) { P mrs (k)P jr (k) - p rs (k)P mr (k)} = 0 
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Thus, this potential divergence of the p integration cancels. The term of first order in q 
vanishes after spherical integration. The second order term will lead to the scaling 

p2 k f pk* p2 k l pk' 

/ dk E(k ) / dq E(q)k(q 2 /k ) ~ / E(k) dk / dq q 2 E(q) ~ (fc') 2/3 
Jk' JO Jk' Jo 

Note that the second term in Eq. (11) also scales this way: the leading order term vanishes 

because it is of odd order and the term of next order scales as above. We conclude that 

J + ~ (k') 2/3 (12) 

Adding Eqs. (9), (10), and (12), 

i ~ -J + + Zf +Z 2 - ~ e 4 / 3 (fe') 2/3 (13) 


So far, the scale k' is arbitrary and i has been evaluated assuming an infinite Kol- 
mogorov inertial range. To apply these results to single point turbulence modeling, it is 
appropriate to replace the infinite Kolmogorov range by a spectrum with a finite Kol- 
mogorov range exending over scales ko < k < kj. The convergence of the integrals under 
consideration implies that they are weak functions of the actual energy spectrum in the 
low and high wavenumber regions k < ko and k > kj. Accordingly, we will make the 
analytically simplest choice, that the spectrum vanishes outside this range, as in Eq. (6). 

To fix the scale k', we reason that in single point turbulence modeling, e is the power 
input into the inertial range, so that i should be evaluated at k' = ko, corresponding to 
the largest inertial range scale. If a model of this type were to be applied in a large eddy 
simulation, k' would be equated instead to the inverse filter size. In the context of single 
point modeling, the choice k' = ko is also consistent with the choice of Kfe as a time 
scale, since e/K e^k 2 ' 3 is the Komogorov inertial range frequency evaluated at ko. 
The turbulence kinetic energy K is given in the infinite Reynolds number limit ko — *• 0 by 

poo o 

K = / dk E(k) = '-C K £ 2/i k~ 2/3 (14) 

Jk 0 2 

Substituting Eqs. (3), (4), and (6) in Eq. (2) and setting k' = ko, the cutoffs cause 
all terms to vanish except for the last term in I ~ , so that 

e = dk Pi mn (k)Pir 3 (k) f dp dq P mr (p)Pns{q)Q(p)Q(q) (15) 

’ Jo </ k=p+q;p , g > fc 0 
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Eq. (13) with k' = ko and Eqs. (14)-(15) then imply that 

£ = -C e 2~ (16) 

Although the constant C e 2 could be evaluated from Eq. (15), our concern here is only to 
derive the forms of the relevant models. Note that the demonstration that the integrals 
converge is necessary to prove the proportionality of i to e 2 /K because dimensional analysis 
alone cannot rule out dependence on Rey. 

III. The dissipation rate equation in rotating turbulence 

Im posed time scale effects like rotation can be included in this analysis provided that 
time scale dependent field descriptors are known. Thus, a rotation dependent equation for 
e will follow by substituting rotation dependent expressions for the correlation Qij in Eq. 
(2). Appropriate expressions appear in the recent work of Mahalov and Zhou. 17 This pro- 
cedure could help address theoretical questions concerning anisotropy of spectral transfer 
in rotating turbulence. In order to evaluate the overall effect of rotation on the dissipation 
rate equation, we will apply Zhou’s simplified isotropic theory of rotating turbulence 15 
based on a phenomenological prescription for the rotation dependent decorrelation time. 
Naturally, geometric effects like two-dimensionalization are not accessible to such a theory; 
a more complete theory would begin with anisotropic expressions and then evaluate shell 
averages. 

The theory of Ref. 15 predicts that for large rotation rates Cl, the energy spectrum is 
given by 

B(p,S2) = C|(e|Q|) 1 /V 2 (17) 

where the modified Kolmogorov constant could be computed from the flux integral for the 
k~ 2 energy spectrum and k° time scale postulated 15 for the strong rotation limit. For 
general rotation rates, 

E(p,Cl) = Z(p,Cl)E(p) (18) 

where E(p) is the rotation independent Kolmogorov spectrum and the rotation dependent 
correction Z is given in Ref. 15. Substituting Eqs. (18), (3), and (4) in Eq. (2) leads to an 
expression for e which cannot be evaluated analytically; the resulting integral could only 
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be given in tabular form as a function of 0. Since the rotation dependent spectrum of Eq. 
(18) is only approximate, it is reasonable to simplify the theory further. 

We will follow a procedure which has now become widespread, 18 namely to evaluate 
i in the strong rotation limit, and then bridge the strong and weak rotation limits by an 
ad hoc interpolation. Substituting the strong rotation spectrum of Eq. (17) in Eq. (2), 

i = - j(Cf) 2 £ \n\j dk Pi„„(k)Pi„(k) 

/ dp dq P‘, mr(jp)P\ rea(q)(p?) ^ (19) 

J k=p+q;p,?>feo 

The integral in Eq. (19) can be shown to be convergent and independent of &o ; therefore 


e = -eg JO ]e 


( 20 ) 


in the strong rotation limit. 

A consequence of this model is the complete suppression of decay of turbulence in 
the strong rotation limit: the solution of Eq. (20) with the energy equation K = —e and 
initial conditions K = Kq , e = £o when t = 0 is 

te " n ‘- 11 

e/e 0 = e~ cm 


where C = C% • The long time limit e = 0 indicates the complete suppression of energy 
transfer in this model. 

A simple formula which interpolates between the strong rotation limit defined by Eq. 
(20) and the case of no rotation defined by Eq. (16) is the Bardina 19 correction 


i = 



C.2 e 1 


But unless helicity is explicitly included in the theory, the corrections for small 0 must 
be independent of the sign of Q and hence quadratic to lowest order. A simple expression 
with this additional property is 



( 21 ) 
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Speziale 16 has recently obtained a model with the same form as Eq. (21) by an entirely 
different argument beginning from stress transport models. Simulations 20 of rotating de- 
caying turbulence using Eq. (21) suggest that fully satisfactory comparison with DNS data 
may require low Reynolds number corrections of the model. 16 

IV. The dissipation rate equation in buoyant turbulence 

In buoyant turbulence, the additional spectral flux e$ appears; the analog of Eq. (1) 
is 

e* « i{I- - I+}k m < T(p)u m (q)T(-k) > (22) 

where T denotes the temperature. The correlations in buoyant turbulence are defined by 

Q t ( k) =< T(k)T(— k) > 

Qr(k) =< T(k)u r (— k) > 

QU k) =< ui(k)uy(-k) > (23) 

Differentiating Eq. (22) with respect to time, following the procedure of Sect. II, 
e e = i{I~ - I + }k m {< T(p)u m (q)T(-k) > + < T(p)u TO (q)T(-k) > 

+ < T{p)um(q)f(-k) >} 

= {/" - I + l[T 1 (k > p ) q) - gT\ k,p,q)] (24) 

where the direct interaction approximation gives 

T 1 = —{k m p n Q T (k)Q mn (q L ) + A; m p n Q^(q)Q^(k) 

+ k m P mra (q)Qr(k)Q h s (p) 

~ kmknQm(P)Qn(q) ~ k m k n Q T (p)Q^ n (q)} (25) 


T 2 = ife m P m 3(q){pr0(p,k,q)[Qr(k)Q :r (q) + Q T (k)Q^(q)] 

+ g r 0(q, k, p) [Qr(k)Q T (p) + Q T (^)Qr (P)] 

- fc r 0(k, p, q) [Q T (p)Qr (q) + Q T {q)Q h r {p)] } (26) 
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and gravity acts in the 3 direction. The constant g denotes the product of the acceleration 
of gravity and the thermal expansion coefficient of the fluid. The T 2 terms originate from 
the u term in Eq. (24); the quantity 0 in Eq. (26) is a wavenumber-dependent time scale. 

Recent theoretical arguments 21 ’ 22 suggest that buoyant turbulence can exhibit Bol- 
giano scaling 23 in which the spectra corresponding to Q T , Q h ,Q u are 

E u = C 9 K g 4/s el /5 &“ u/ ® 

E h = Clg 1 ** 4 /5 fe- 9/5 

E t - dg- 2/5 ey 5 k~ 7 / 5 (27) 

In this scaling, the time scale © has the form 

@(k, P, q) = | n(k) + ar/(p ) + rj(q)] _1 (28) 

where 

<m = cvV* 2 ' 5 ' (29) 

and a is the inverse turbulent Prandtl number for free convection. A preliminary attempt 
to compute the inertial range constants a,Cj <: ,Cl,Cj’, and Cf, appears in Ref. 24. Bolgiano 
scaling should describe buoyant flows, or at least scales of motion, in which the gravitational 
coupling dominates nonlinearity, although the occurence of this scaling in Rayleigh-Benard 
convection remains controversial. 25 If instead, the gravitational coupling is small compared 
to the nonlinearity, then the temperature is a passive scalar with a source of fluctuations; 
Kolmogorov scaling applies, in which 

E u = C K s 2/3 k~ 5/3 
E h = 0 

E t = C B £9£~ 1/3 k- 5 / z (30) 

Kolmogorov scaling applies to locally isotropic turbulence; accordingly, the heat transfer 
spectrum E h vanishes in Eq. (30). Bolgiano scaling applies to anisotropic turbulence in 
which the direction of gravity is distinguished; this permits a nonvanishing E h in Eq. (27). 

The dissipation rate equation for each type of scaling will be treated separately. 
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A. Bolgia.no scaling 

In evaluating the integral Eq. (24) with Eqs. (25)-(27), the isotropic contribution to 
the energy transfer will be evaluated by setting 

Q t ( k) = Q T (k) 

Qr( k) = Q h (k)P r 3 (k) 

er,(k) = o*(*)«j(k) 

When these substitutions are made, the index 3 always appears in pairs. Thus, to evaluate 
the isotropic part of the energy transfer, it is convenient to sum over this index. This 
procedure leads, as in Sects. II and III, to an integral for the proof that this integral 
converges is similar to the proof in Sect. II and therefore need not be given. The terms 
T 1 scale after integration as 

{/" - I + }T 1 ~ g 2/5 e 6 9 /5 (k') 2/5 (31) 

Integrating E T from Eq. (27) over a Bolgiano spectrum with lower cutoff ko leads to 

(32) 

thus, in terms of single point quantities, setting k' = ko 

{I~ - I+JT 1 = -J f° dk k m k n f dp dqx 

^ JO </k=p*fq;p,g>^o 

{5<3S.,(p)0^(q) + Q T (p)Qln(<i)} (33) 

or in view of Eqs. (31)-(32), 

{/+ - r }T' = (34) 

Analogous treatment of the T 2 terms gives 

{J--/ + }T 2 ~ 9 - 3/5 4 /5 fco /5 (35) 

Integrating the Bolgiano spectra E h and E u of Eq. (27), 

h = f c'/ /5 4 /! V /5 

K = lc e K g >/S 4 ,5 ^ /b (36) 
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where the heat transfer H is the single point moment H = — il^T >. Substituting Eq. 
(36) in Eq. (35), 

{/- - I + }T 2 = Cl^H (37) 

Combining the results of Eqs. (34) and (37), 

it = - CUaj: < «*T > (38) 

which agrees with the model proposed in Ref. 2 from which the nomenclature for the 
model constants is taken. 

B. Kolmogorov scaling 

Substituting the Kolmogorov spectra Eq. (30) in Eqs. (24)-(26) leads to 

{/--J+JT 1 ~ — ege 

The term T 2 vanishes since Q h = 0; therefore, 

ie = -Cl^e e (39) 

Comparison of Eqs. (38) and (39) reveals that Bolgiano and Kolmogorov scaling 
have fundamentally different implications for modeling buoyant turbulence. First, the g 
dependent term of Eq. (38) is absent in Eq. (39). More fundamentally, the time scale in 
Eq. (38) is the thermal time scale Kq/sq whereas the time scale in Eq. (39) is the velocity 
time scale K/e. The results for different scaling laws are combined in the model of Ref. 2, 
in which the destruction terms are 

id = - CLj^ed - < u,T > (40) 

A model of this form is not inconsistent with the present theory, since Bolgiano and Kol- 
mogorov scaling can both occur in the same problem at different inertial range scales. 21 ’ 22 
However, particular problems can be dominated by only one type of scaling: computational 
experience indicates 26 that the g dependent term in Eq. (40) improves computational 
predictions for some problems but degrades them in others; this difference may reflect 
dominance in these problems by different scaling laws. Future work must clarify how the 
existence of these two scaling laws should be reflected in turbulence models. 
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V. Rotating buoyant turbulence 

A modeled transport equation for eg in rotating buoyant turbulence follows from 
Eqs. (24)-(26) by substituting rotation-dependent field quantities Q r , Q\ and Q?-; the 
calculation does not differ in principle from that of Sect. Ill for rotating incompressible 
turbulence. As in Sects. Ill and IV, the present analysis will be limited to a simplified 
isotropic analysis which treats only the overall effect of this type of coupling; the subtle 
geometric effects which can arise if the rotation axis does not coincide with the direction of 
gravity are inaccessible to this type of theory. The discussion will be limited to Bolgiano 
scaling. 

The approximate analysis of Bolgiano scaling in Ref. 24 expresses E h and E u in terms 
of E t by 

E h ~ gQE T E u ~ (gQ) 2 E T (41) 

where 0 is a time scale. A rotation dependent spectrum E T can be found from the DIA 
flux balance for Bolgiano scaling. For the purpose of computing scalings, this balance takes 
the form 

eg ~ g 2 k A QE T E u (42) 

Following Zhou 15 , assume that 0 ~| 0 | -1 in the limit of strong rotation. Substituting E u 
from Eq. (41) in Eq. (42), 

E t = C%°g~ l | Q | 3/2 e x J 2 k ~ 2 (43) 

Then Eq. (41) implies 

E h = C% e j 0 | 1/2 e l J 2 k~ 2 

E u = C^g | Q |- 1/2 ej 2 k~ 2 (44) 

It is crucial that the DIA flux integral, given in simplified form as Eq. (42), converges when 
the spectra of Eqs. (43)-(44) are substituted in it. The convergence of the flux integral is 
easily checked; only because of this convergence can the scaling law of Eq. (43) be inferred 
from Eq, (42). 

Substitution of Eqs. (42) and (44) in Eq. (24) leads as in Sect. Ill to the strong 
rotation limit 

eg = -C™\n\eg (45) 
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An ad hoc bridge between the weak and strong rotation limits can be constructed as in 
Sect. III. A possible form is 


6$ == _ Ce2~i^~ £e {]- ~*~ 


CgilKi^ 1/2 

0% U J ' 


- C! f gjt < u 3 T > {1 + 


eg ftg .2,1/2 

d.geeH 1 ’ 


VI. Conclusions 

We s umm arize the models obtained for the destruction terms in dissipation rate trans- 
port equations: 

A. Incompressible turbulence 

t-c.4 

where C e 2 can be computed from Eq. (15). Convergence of this integral at high wavenum- 
bers implies that C e 2 does not depend on turbulent Reynolds number Rer- 

B. Rotating turbulence 

For rapid rotation, 

i = -C% | ft | e 

where C% is found from Eq. (19). A possible interpolation formula between the weak and 
strong rotation limits is 

s = -C e2 — { 1+ [ 

C. Buoyant turbulence 

For Bolgiano scaling, 

i» = - C! 5 gjt < u 3 T > 

For Kolmogorov scaling 

The time scale for Bolgiano scaling is ee[K$\ the time scale for Kolmogorov scaling is e/K. 
The buoyancy effect present for Bolgiano scaling is absent for Kolmogorov scaling. 


C% fog -. 2 -j 1/2 
C S 2 e n 
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Bolgiano and Kolmogorov scaling are believed to occur together in general flows; this 
suggests that the general dissipation rate model could combine the results for both scalings 
as in Eq. (40). 2 However, a complete model must provide for flows dominated by either 
scaling law alone. 

D. Rotating buoyant turbulence 

In the strong rotation limit, for Bolgiano scaling 

io = -C% \a\e. 

As in the case of rotating incompressible turbulence, an ad, hoc expression interpolating 
the weak and strong rotation limits can be constucted. 
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